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摘 要 : 厘清 气候 变化 和 人 类 活动 对 植被 净 初 级 生产 力 (NPP) 的 影响 ,对 区 域 生态 系统 变化 与 可 持续 发 展 具 有 重 
要 意义 。 基 于 CASA 模型 ,利用 1994 一 2018 年 5 期 Landsat 遥 感 影像 及 气象 数据 ,通过 仿真 情景 实验 研究 了 人 类 活 
动 和 气候 变化 对 哈萨克 斯 坦 首都 圈 植被 净 生 产 力 的 影响 。 结 果 表 明 :(1) 1994 一 2018 年 间 , 哈 萨克斯 坦 首都 圈 植 被 
NPP 是 波动 上 升 趋势 ,NPP 多 年 均值 为 226.21 g Cm*…a';(2) 首都 圈 人 工 造林 面积 不 断 增加 ,对 NPP 有 显著 的 正 增 
益 效 果 (0.38 Tg Cea! P<0.01) ,而 气候 变化 对 NPP 影 响 波 动 比较 大 ,整体 上 呈现 负 增 益 效果 (-0.07 Tg C+ a’, P= 
0.34) ;在 土地 履 盖 变化 和 气候 变化 共同 作用 下 ,哈萨克 斯 坦 首都 圈 植 被 NPP 总 体 上 呈 显 著 正 增益 效果 (0.27 Tg Cra, 
P<0.1)。 相 较 于 气候 变化 ,人 类 活动 对 植被 净 初 级 生产 力 影响 更 加 明显 ,相关 性 也 更 强 ;(3) 气候 因子 中 温度 ,太阳 
辐射 和 降水 是 影响 NPP 的 主要 因素 。1994 一 2000 年 和 2006 一 2012 年 期 间 由 于 温度 的 升 高 和 降雨 量 的 减少 导致 
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NPP 损 失 , 分 别 日 


H 218.50 g C+m*+a' .201.19g Cm a 


降 到 189.00 g C-m*+a'.188.48 g C -m° a” ;2000 一 2006 年 和 


2012 一 2018 年 期 间 随 着 降水 条 件 的 改善 ,该 地 区 NPP 均 值 显著 上 升 ,均值 分 别 达 到 201.19 g C+m?+a' 207.73 g C- 
ma"'。 本 文 研究 结果 有 助 于 揭示 气候 变化 和 人 类 活动 对 植被 NPP 的 影响 机 制 ,可 为 改善 荒漠 草原 区 生态 环境 质 


中 | 


量 , 缓 解 全 球 变 暖 问题 和 实现 哈萨克 


斯 坦 2060 年 碳 中 和 目标 提供 参考 。 


Kia: CASA 模 型 ，NPP; 哈萨克 斯 坦 首都 轿 ; PCA 分 析 ; 仿真 实验 设计 


植被 净 初 级 生产 力 (Net Primary Productivity , 
NPP) 是 指 植物 在 单位 时 间 .单位 面积 内 所 积累 的 有 
机 物 的 量 ,是 由 植物 光合 作用 所 产生 的 有 机 质 总 量 
(Gross Primary Productivity ,GPP) 中 减 去 自 养 呼吸 后 
的 剩余 部 分 ,也 称 为 净 第 一 生产 力 "。NPP 被 认为 
是 生态 系统 过 程 的 调节 需 和 陆地 碳 汇 的 主要 决定 
因素 ”, 而 陆地 生态 系统 是 目前 全 球 碳 循环 的 主要 
汇 , 湾 在 地 抵消 大 量 人 为 碳 排放 。 因 此 ,NPP 不 仅 
直接 反映 植物 群落 在 自然 环境 条 件 下 的 生产 能 力 ， 
表征 陆地 生态 系统 的 质量 状况 ,而 且 是 测算 土地 履 
盖 及 能 源 消 费 碳 排放 ,判定 生态 系统 碳 源 . 汇 以 及 
调节 生态 过 程 的 主要 因子 “ ,在 全 球 变化 及 碳 平 衡 
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中 起 着 至 关 重 要 的 作用 ”” 。 

由 于 地 面 测量 数据 无 法 描述 NPP 在 较 大 区 域 
范围 内 的 时 空 分 布 特征 ,所 以 使 用 多 源 遥 感 数据 和 
数学 模型 来 模拟 计算 NPP 就 变 成 了 一 种 很 常用 的 
重要 方法 。 相 比 于 传统 的 生态 学 过 程 模 型 和 气候 
生产 力 模 型 涉及 参数 多 旦 复杂 度 高 和 偏差 较 大 的 
缺点 , 光 能 利用 率 模型 数据 获取 难度 较 小 同时 模型 
模拟 精度 较 高 ,还 可 以 通过 遥感 手段 辅助 计算 ,节省 
了 大 量 繁琐 的 野外 试验 测定 步骤 ,成 为 了 NPP 佑 算 
中 的 主要 研究 方法 “” 。 其 中 ,CASA 模 型 利用 归 一 
化 植被 指数 (NDVI) 和 土地 覆盖 分 类 2 种 遥感 数据 ， 
以 及 结合 月 平均 温度 、 太 阳 辐 射 .降水 SEFA 
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象 数据 对 NPP 进 行 估算 ,在 减少 复杂 参数 收集 的 同 
时 也 可 以 更 加 精确 地 完成 大 范围 的 NPP 模 拟 7” ， 
在 国内 外 被 广泛 应 用 于 不 同 区 域 各 种 地 表 植 被 类 
型 固 碳 量 的 反 演 模拟 一 。 

哈萨克 斯 坦 首 都 阿 斯 塔 纳 地 处 中 亚 干旱 - 半 王 
时 气候 带 '“”, 是 世界 上 第 二 冷 的 首都 城 市 , 且 大 
风 沙尘暴 及 暴风 雪 频 发 ,生态 系统 极为 脆弱 2 。 
因此 ,研究 哈萨克 斯 坦 首 都 圈 植 被 NPP 时 空 变化 ， 
探究 气候 变化 与 人 类 活动 影响 下 的 驱动 因素 ,厘定 
气候 因素 中 的 主导 因子 ,评估 气候 变化 与 人 类 活动 
对 亚 寒带 草原 生态 系统 NPP 的 影响 ,对 于 区 域 生态 
系统 的 可 持续 管理 起 着 至 关 重 要 的 作用 。 


1 研究 区 概况 


哈萨克 斯 坦 首都 阿 斯 塔 纳 ( 努 尔 苏 丹 ,51.00?~ 
51.30°N ,71.22°~71.74°E, 图 1) 位 于 辽阔 的 哈 陛 克 大 
草原 中 北部 ,属于 典型 温带 大 陆 性 气候 区 。 年 平均 
气温 3 % ,平均 海拔 374 m, 多 年 平均 降水 量 290 
mm。 冬 季 寒 冷 漫长 ,暴风 雪 严 重 ,夏季 短 热 淹 湿 。 
自 1997 年 迁都 至 今 ,首都 圈 人 口 .经济 以 及 城区 建 
设 得 到 了 飞速 发 展 ,如 今 阿 斯 塔 纳 已 经 成 为 仅 次 于 
前 首都 阿拉 木 图 的 哈萨克 斯 坦 第 二 大 城市 。 
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2 数据 与 方法 


2.1 数据 来 源 

ASCH GEE (Google Earth Engine) 平 台 获 取 CA- 
SA 模型 所 用 到 的 有 归 一 化 植被 指 数 (Normalized 
Difference Vegetation Index, NDVI) .土地 和 覆盖 类 型 
数据 和 温度 、 太 阳 辐 射 量 、 实 际 菩 散 量 潜在 蒸 散 量 
等 I。 其 中 ,NDVI 是 对 Landsat 数 据 再 次 去 云 处 理 
并 进行 波段 运算 得 到 ,空间 分 辩 率 为 30 m。 土 地 覆 
盖 类 型 数据 是 基于 1994 一 2018 年 的 Landsat 影像 ， 
参考 全 球 30 m 土 地 覆盖 数据 (GlobeLand 30) 的 分 类 
标准 及 制图 规范 ,建立 了 哈萨克 斯 坦 首都 圈 土 地 禾 
盖 分 类 系统 。 在 此 基础 上 ,利用 随机 森林、 目 视 解 
译 相 结合 的 方法 ,对 哈 了 村 克 斯 坦 首 都 阿 斯 塔 纳 周 边 
区 域 的 土地 覆盖 类 型 进行 遥感 分 类 ,经 过 严格 的 修 
改 和 精度 验证 ,得 到 较为 精确 的 土地 覆盖 分 类 

Terra Climate 数据 来 自爱 达 答 大 学 发 布 的 全 球 
陆地 表面 月 均 气 候 和 气候 水 平衡 数据 集 ,涵盖 了 全 
球 气候 水 平衡 和 陆地 表面 每 月 气候 数据 。 包 括 太 
阳 辐 射 .潜在 蒸 散 量 .实际 蒸 散 量 等 气温 数据 。 采 
用 MODIS 数 据 集 中 的 月 平均 温度 数据 以 及 NPP 产 
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图 例 
。 阿 斯 塔 纳 ( 努 尔 苏 丹 ) 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 CS(2021)5453 号 ,对 底 图 边界 无 修改 。 下 同 。 


图 1 阿 斯 塔 纳 及 周边 地 区 主要 土地 覆盖 类 型 


Fig. 1 Main land cover types in and around Astana 
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品 , 其 中 NPP 产 品 用 于 对 CASA 模型 模拟 结果 进行 
精度 验证 。 利 用 GEE 云 平 台 下 载 上 述 气象 数据 进 
行 重 投影 裁剪 操作 ,用 于 CASA 模型 的 计算 。 
2.2 研究 方法 
TE CASA 模型 中 ,植被 吸收 的 有 效 光 合 辐射 以 
及 在 进行 光合 作用 过 程 中 实际 的 光 能 利用 效率 直 
接 决 定 了 该 模型 对 NPP 的 模拟 计算 结果 。 本 文 利 
用 朱文 泉 等 * 在 干 量 区 生态 研究 中 已 经 得 到 应 用 
和 验证 的 研究 成 果 , 配 置 CASA 模型 最 大 光 能 利用 
CASA 模 型 中 的 NPP 主 要 由 植物 的 光合 有 效 辐 
射 (APAR) 和 实际 光 能 利用 率 (s)2 个 因子 来 控制 ， 
NPP(x,t) = APAR(x,t) x e(x,t) (1) 
植被 光合 作用 过 程 中 吸收 的 有 效 光 和 辐射 
(APAR ) 取 决 于 太阳 总 辐射 量 以 及 植被 对 有 效 光 合 
辐射 的 吸收 比例 ,计算 公式 如 下 
APAR(x,t) =SOL(x,1)x FPAR(x,t)x0.5 (2) 
式 中 : APAR(x 沁 表示 某 像 元 x 在 上 月 吸收 的 光合 有 
效 辐射 (g Cem HT); ex 为 表 示 某 像 元 x 在 上 月 的 
实际 光 能 利用 率 (g CMI’) , 本文 参考 朱文 泉 等 " 
研究 成 果 ,确定 了 不 同 植被 类 型 最 大 光 能 利用 率 
( 表 1)。 SOLS, ) 表示 1 月 在 像 元 x 处 的 太阳 总 辐 
射 量 (g Cem 4!) ; FPAR(«,0) 表示 植被 层 对 于 入 
射 有 效 光 合 辐射 的 吸收 比例 ;常数 0.5 表示 植被 进 
行 光 合作 用 过 程 中 所 能 利用 的 太阳 光 有 效 辐射 波 
段 范围 (波长 为 0.4~0.7 um ) 占 太阳 总 辐射 量 的 比 
例 。 


表 1 不 同 植被 类 型 最 大 光 能 利用 率 
Tab. 1 Maximum light energy utilization rate of different 


vegetation types 


编号 植被 类 型 < 模拟 值 
1 落叶 针 叶 林 0.485 
2 ER ETAT PK 0.389 
3 TEI Ra} N 0.692 
4 TER AAP PK 0.985 
5 针 间 混交 林 0.475 
6 HER TAT TR EE SCA 0.768 
7 HEMA 0.429 
8 草地 0.542 
9 耕地 0.542 
10 其 他 0.542 
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环境 因子 如 大 气温 度 .土壤 含水 量 等 通过 影响 
植物 的 光合 能 力 , 即 光 能 利用 率 e(x,t) 来 调节 植被 
的 NPP。 


&(x,t) = T aien x T 92.) x W test) XE a (3) 

2 
T „a =0.8 +0.02 XT, oy ~ 0.0005 x [Pano] (4) 

_ 1.184 y 

€2(x,1) 
{1 + exp[0.2%7,,4.)-10- 7.0] 
(5) 
1 
{I + exp[0.3 x CT, y)— 10- Tal 

Wo.) =0.5 +0.5 x EET, /EPT (6) 


式 中 : Taea 和 7 表示 低温 情况 和 高 温情 况 对 光 
能 利用 率 的 胁迫 作用 ; Wy, 表示 水 分 胁迫 影响 系 
数 , 反 映 水 分 条 件 的 影响 ; s,, 表示 理想 条 件 下 的 
最 大 光 能 利用 率 (g CMI); Ta 表示 植物 生长 的 
最 适 温度 ; Ty) 表示 月 平均 温度 (% );EET 表 示 区 域 
实际 蒸 散 量 (mm);EPT 为 区 域 潜 在 蒸 散 量 (mm)29。 
2.3 仿真 实验 设计 

本 文 设计 了 3 种 仿真 情景 来 模拟 1994 一 2018 
年 气候 变化 及 人 类 活动 引起 的 土地 覆盖 类 型 变化 ， 
评估 其 对 哈萨克 斯 坦 首都 圈 植 被 NPP 的 影响 ( 表 
2)。 人 情景 1: 按 照 土 地 覆盖 类 型 和 气候 变化 的 实际 
情况 ,利用 相应 年 份 的 土地 覆盖 类 型 和 气候 数据 对 
植被 NPP 进 行 模拟 计算 ;情景 2: 控 制 土地 覆盖 类 型 
和 NDVI 保 持 在 1994 年 的 水 平 ,研究 从 1994—2018 
年 以 来 气候 变化 对 植被 NPP 的 影响 ;情景 3: 控 制 气 
候 数 据 保 持 在 1994 年 的 水 平 ,研究 1994 一 2018 年 
人 类 活动 土地 覆盖 变化 对 植被 NPP 的 影响 。3 种 情 
景 下 对 NPP 的 模拟 均 利 用 Python 语 言 计 算得 到 。 

为 了 进一步 探究 温度 降水 太阳 辐射 AE: 
各 气候 因子 对 植被 NPP 的 影响 ,随机 选取 500 个 点 


表 2 不 同 模拟 仿真 情景 


Tab.2 Different simulation scenarios 


仿真 情景 LULC NDVI sol temp aet pet 
土地 覆盖 类 型 和 气候 变化 
土地 芥 羡 类 型 变化 A A O O oo O 
气候 变化 Oo Oo 


注 :人 表示 输入 量 从 1994 一 2018 年 的 变化 ;O 表 示 输 入 量 保持 不 
变 ;LULC 表 示 土 地 覆盖 类 型 ;NDVI 表 示 归 一 化 植被 指数 ;sol 表 示 
太阳 辐射 ;temp 表示 平均 温度 ;aet RIR SE PNAS BG ; pet FAN TATE ZR 
散 量 。 
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进行 主 成 分 分 析 和 皮尔 森 相关 系数 计算 ,进一步 研 。” 合理 的 且 适 用 于 此 区 域 。 

究 NPP 变 化 的 气候 驱动 因子 。 

2.4 模型 评价 3 结果 与 分 析 


空间 异 质 性 .时间 上 的 不 一 致 性 以 及 缺乏 直接 
测量 的 NPP 值 ,如 何在 大 范围 内 建 模 验证 NPP 一 直 
是 陆地 生态 碳 积累 研究 中 的 重点 … 。MODIS 在 进 
行 NPP 计 算 的 时 候 对 很 多 区 域 进行 了 插值 补充 ,在 
研究 大 范围 NPP 的 情况 下 该 产品 具有 较 好 的 应 用 ， 
但 是 无 法 满足 在 更 精细 尺度 下 探讨 NPP 的 驱动 因 
素 。MOD17A3H V6 产品 提供 了 涵盖 整个 中 亚 地 区 
500 m 空间 分 辨 紊 的 年 度 净 初 级 生产 力 (NPP)。 
此 ,本 文 利用 MODIS 的 NPP 产 品 来 验证 CASA 模型 
模拟 整个 地 区 NPP 结果 的 合理 性 。 

验证 结果 表明 :CASA 模型 模拟 哈萨克 斯 坦 首 
都 圈 不 同 土地 覆盖 类 型 的 平均 NPP 4 MODIS 数据 
产品 一 致 性 很 好 ,特别 是 防护 林 和 裸 地 (图 2a) ,在 
整个 哈萨克 斯 坦 首都 圈 MODIS-NPP 与 CASA-NPP 
的 平均 值 均 具 有 较 好 的 相似 度 (P<0.001)。 模 型 模 
拟 结果 显示 (图 2b), 在 种 植 防护 林 的 区 域 NPP 值 明 
显 增加 ,最 大 变化 值 达到 318.7 g Com”, Hhh NPP LY 
值 变 化 不 大 , 裸 地 和 水 体 以 及 建筑 区 NPP 值 出 现 明 
显 下 降 , 最 大 下 降 300.9 gC.m*( 图 2b)。 总 体 分 析 
来 看 ,该 地 区 NPP 保 持 不 变 的 区 域 占 比较 大 。 基 于 
以 上 验证 结果 ,本 文 认 为 CASA 模型 模拟 的 结果 是 


(a) 


MODIS-NPP/(g Cm?:a-!) 


CASA-NPP/(g C-m*-a"') 


3.1 蛤 萨克斯 坦 首都 圈 NPP 变化 特征 

基于 CASA 模 型 模拟 的 结果 ( 表 3), 随 着 研究 区 
土地 上 覆盖 类 型 的 改变 ,防护 林 面 积 不 断 增加 (图 3a， 
图 3b) , 1994—2018 年 首都 轿 植 被 NPP 呈现 先 下 降 
后 增长 的 波动 上 升 趋势 ,NPP 年 均值 为 226.21 gC- 
m2。 其 中 ,1994 一 2000 年 NPP 均 值 每 年 下 降 约 
2.74 g C.m2.a ,从 218.50gC.om2 下 降 至 202.07 g 
Cm”, 随 后 则 逐年 上 升 ;2006 年 .2012 年 NPP 均值 
分 别 为 216.72 g Com? 45 235.35 g C+m”, 2018 42 NPP 
均值 为 258.42 g C.m’, 达 到 该 地 区 最 大 值 。 在 
1994—2018 年 间 , 首 都 圈 内 城镇 建筑 水 体 和 裸 地 
区 域 NPP 出 现 明显 的 损失 (图 3c) ,最 高 损失 达到 
325.28 g Cm”, 在 其 他 植被 和 人 工 林 区 域内 NPP 均 
值 显著 增加 ,特别 是 人 工 林 区 域 NPP 平 均值 增加 最 
为 显著 ,最 高 达到 310.18 g Cem. 

人 研究 发 现 , 在 1994 一 2018 年 期 间 , 首 都 圈 防 护 
林 、 草 地 、 裸 地 、 建 筑 用 地 及 水 体 等 不 同 土地 覆盖 类 
型 下 植被 NPP 均 值 由 大 到 小 为 :防护 林 > 其 他 植被 > 
裸 地 > 建筑 用 地 > 水 体 。5 种 土地 覆盖 类 型 NPP 年 均 
值 从 大 到 小 为 :防护 林 > 草 地 > 裸 地 > 建筑 用 地 > 水 


(b) 


人 


NPP 变 化 值 
Ag Cm?) 


= . 
0 


e—a 


NPP/(g Cm?) 


图 2 CASA 模 型 精度 检验 与 研究 区 NPP 变 化 
Fig.2 CASA model simulation NPP accuracy verification 
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表 3 各 土地 覆盖 类 型 不 同年 份 NPP 值 
Tab.3 NPP values of different land cover types in different years 


(a) 1994 年 NPP 


>z 


561.25 


(c) 1994 一 2018 年 NPP 变 化 


0 15km 
[= 


NPP/(g C-m2) 


—300.94 


年 份 防护 林 /(g Cem?) 草地 /(g Cem?) 水 体 /(g Com?) 建筑 区 /(g Com) 裸 地 /(g Com?) 地 区 年 均值 /总 值 /(g Cm’) 
1994 均值 317.11 195.73 150.75 185.65 198.34 218.50 
总 值 3.78x 10" 1.93x 10° 5.81x 10’ 1.16x 10’ 1.53x 10° 3.56% 10° 
2000 均值 272.27 186.61 187.09 198.66 167.06 202.07 
总 值 6.81x 10" 2.00% 10° 5.49x 10’ 2.82x 107 1.12x 10° 3.27% 10° 
2006 均值 214.10 198.76 125.12 204.67 192.68 216.72 
总 值 1.15x 10° 1.63x 10° 3.74x 10” 2.77 10’ 1.72x 10° 3.53x 10° 
2012 均值 266.13 212.43 191.42 196.88 206.60 235.35 
总 值 2.31x 10° 1.74x 10° 4,32x 10’ 2.86x 10’ 1.80% 10° 3.83% 10° 
2018 均值 302.81 241.14 189.98 163.86 216.99 258.42 
总 值 2.94x 10° 2.17 10° 7.53x 10" 2.32 10" 1.65x 10° 4.21x 10° 


(b) 20184 NPP 


>z 


NPP/(g C-m?) 
E 580.87 
37.25 
9 isim 
(d) NPP 年 均值 变化 
— 实际 值 
280 o 气候 变化 


一 - 人 类 活动 变化 
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NPP/(g Cm >a) 
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图 3 人 类 活动 与 气候 共同 作用 下 的 首都 圈 NPP 结 果 


Fig. 3 NPP results under the joint action of human activities and climate 


体 , 相 应 地 5 种 土地 覆盖 类 型 占 地 面积 分 别 从 1994 
年 的 1.02% 、58.47% 、37.87% .0.35% 、2.29% 变 化 至 
2018 年 的 5.35% 、49.68% 、42.00% 、0.78% 、2.19%。 
包括 农田 和 草原 的 其 他 植被 土地 覆盖 类 型 占 比 最 
大 的 , 占 地 面积 由 58.47% 下 降 到 49.68% 植 被 NPP 平 
均值 介 于 186.61~541.14 g Cm 之 间 。 首 都 圈 防 护 
林 面 积 显著 增加 ,防护 林 NPP 均 值 也 达到 最 高 , 介 


于 214.10~317.11 g Cm”, 

由 图 3d 可 知 ,1994 一 2000 年 间 研 究 区 NPP 出 现 
了 损失 ,首都 圈 人 工 造林 的 同时 出 现 NPP 均 值 下 降 
与 预期 目标 不 符合 。 结 合 植 被 NPP 结果 和 土地 利 
用 覆盖 类 型 , 究 其 原因 发 现 ,1994 一 2000 年 时 段 内 
水 体面 积 急 剧 减 少 ,仅仅 6a 间 减少 面积 达 8.24x10" 
m”。 这 也 正 是 自 1997 年 哈萨克 斯 坦 迁都 到 阿 斯 塔 
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纳 后 ,城市 化 建设 导致 短 时 间 内 用 水 过 度 ,水 体 的 
减少 直接 或 间接 导致 了 整个 地 区 NPP 的 减少 。 


2000 年 后 NPP 又 逐渐 上 升 , 很 大 程度 上 依赖 于 首都 
圈 绿 环 建设 工程 的 实施 ,整体 提升 了 NPP 变化 
幅度 。 
3.2 哈萨克 斯 坦 首都 圈 NPP 的 影响 因素 

综 上 所 述 ,25 a 间 首 都 圈 NPP 平 均值 增加 了 
39.92 g Cm?”。3 种 仿真 实验 情景 下 的 NPP 模拟 结 
果 对 比如 图 3d 所 示 , 人 类 活动 对 NPP 有 非常 明显 的 
增益 效果 ,1994 一 2018 年 间 研 究 区 内 NPP 均值 增加 
T 56.45 g Cm ,而 气候 变化 对 NPP 的 影响 呈现 出 
很 大 的 波动 性 ,每 个 年 份 有 增 有 减 ,与 气候 变化 尤 
其 是 降水 变化 趋势 相 吻合 (P=0.618)。 在 人 类 活动 
和 气候 变化 共同 作用 下 ,哈萨克 斯 坦 首都 圈 绿 环 工 
程 的 实施 ,使 得 造林 面积 不 断 增 加 ,对 NPP 有 显著 
的 正 增益 效果 [0.27 Tg Cea (P<0.1)]. 
3.2.1 人 类 活动 变化 对 首都 圈 NPP 影 响 通过 对 比 
图 4a .图 4b 显示 ,1994 一 2018 年 防护 林 面 积 增加 了 
4.68%( 图 4c) ,城市 建筑 区 面积 增加 ,大 部 分 是 通过 
扩张 裸 地 转变 而 来 。 城 区 建筑 面积 所 占 比例 较 小 


(a) 1994 年 土地 覆盖 


>z 


mm 草地 (58.47%) 
ma 水 体 (2.29%) 
æ 建筑 区 (0.35%) 
mm 裸 地 (37.87%) 
0 15km 


(c) 2018 年 人 工 林 空间 分 布 
Fi 


>z 


me 绿 环 工 程 (1.02%) 


各 土地 类 型 面积 占 比 /% 
S 


-a 城市 扩张 对 地 区 NPP 的 影响 极其 有 
民 。 其 他 植被 面积 有 较为 明显 的 下 降 , 绿 环 工程 的 
a 的 灌木 从 /草地 转化 为 森林 (图 4d)。 
防护 林 面 积 的 增加 具有 防风 固沙 的 效果 ,同时 
能 够 提高 地 区 陆地 生态 系统 固 碳 能 力 以 及 改善 当 
地 生态 环境 质量 。 基 于 情景 2, 人 类 活动 (土地 覆盖 
类 型 变化 ) 对 首都 图 陆地 净 生 产 力 具有 积 ! 极 的 影响 
(图 $) , 主要 得 益 于 和 森林 面积 的 增加 (图 4c)。 在 该 
情景 下 首都 圈 陆 地 生态 系统 年 净 生 产 力 增加 显著 
(0.38 Tg C-a', P<0.1) , NPP EMAM 3.56 Tg C 增 加 到 
4.49 Tg C。 其 中 防护 林 、 其 他 植被 水体、 建筑 用 地 、 
裸 地 分 别 贡 献 了 0.26 Tg C、0.33 Tg C、0.02 Tg C, 
0.01 Tg C,0.31 Tg Co 
3.2.2 气候 变化 对 首都 圈 NPP 影 响 1994—2018 4F 
间 首 都 圈 气候 呈现 较 大 的 波动 性 ,但 整体 变化 趋势 
不 明显 (图 6)。 太 阳 辐 射 波动 幅度 较 大 旦 有 上 略微 下 
降 ( 图 6b), 为 -0.019 Wm*…a'; 年 平均 降雨 量变 化 
相对 平缓 ,为 0.14 mm:m?…a'( 图 6a);2000 年 平均 
Bale 25 a 间 波动 幅度 较 
大 ,总 体 上 升 了 0.24 Cea", 


(b) 2018 年 土地 覆盖 


>z 


me 绿 环 工程 (5.35%) 
mm 草地 (49.68%) 
me 水 体 (2.19%) 

mm 建筑 区 (0.78%) 
mm 裸 地 (42.00%) 
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(d) 各 年 份 土地 覆盖 占 比 
we 防护 林 æ 水 体 
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图 4 研究 区 土地 覆盖 类 型 概况 
Fig. 4 The summary of the land cover types within the research region 
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(a) 1994 年 NPP (b) 2018 年 NPP 
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(c) 1994 一 2018 年 NPP 变 化 
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图 5 人 类 活动 作用 下 首都 轿 NPP 变 化 
Fig. 5 Changes in the metropolitan region’ s NPP as a result of human activities 
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(d) 人 类 活动 作用 下 NPP 年 际 变化 
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图 6 研究 区 植被 NPP 与 不 同 气候 因子 变化 情况 
Fig. 6 NPP of vegetation and changes of different climatic factors in the study area 
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气候 因子 如 气温 .土壤 水 分 状况 以 及 大 气 水 汽 
压 差 等 会 通过 影响 植物 的 光合 能 力 而 调节 植被 的 
NPP。 基 于 情景 3 研究 发 现 (图 7) ,气候 变化 对 首都 
圈 NPP 影 响 波动 非常 大 ,但 整体 上 气候 因素 对 NPP 
影响 并 不 明显 ,NPP 均 值 变 化 较 小 (图 7c)。 研 究 区 
西南 部 和 东北 部 NPP 出 现 较为 明显 的 损失 ,西北 部 
NPP 值 有 显著 增加 ,大 部 分 均 保 持 较 微 弱 的 变化 。 

在 温度 胁迫 和 水 分 胁迫 作用 下 ,低温 和 高 温 时 
高 呼吸 消耗 将 降低 光 能 利用 率 , 除 此 之 外 环境 中 植 
物 所 能 利用 的 有 效 水 分 的 增加 也 会 使 光 能 利用 率 
增加 。1994 一 2000 年 和 2006 一 2012 年 间 温 度 的 升 
高 ,尤其 是 降雨 量 的 减少 是 地 表 水 分 减少 的 主要 因 
素 ,直接 导致 了 NPP 损 失 ( 图 7d) ,NPP 均 值 分 别 由 
218.50 g C+m?+a!',201.19 g C+m™°- a” KEZIT 189.00 
g C+m?+a'. 188 g C-m?+a', 2000—2006 年 和 
2012—2018 年 期 间 随 着 降水 条 件 的 改善 ,该 地 区 
NPP 均 值 显著 上 升 ,均值 分 别 达 到 201.19 g Cem? a", 
207.73 g Cm ”a ,气候 变化 对 该 地 区 NPP 整 体 上 
呈现 负 增 益 效果 (-0.07Tg Cra, P=0.34)。 
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(c) 1994 一 2018 年 NPP 变 化 
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3.3 影响 NPP 变化 的 气候 主 控 因 子 

为 进一步 探究 不 同 气候 因子 对 NPP 的 影响 ,在 
该 地 区 影像 上 随机 采取 500 个 点 并 提取 NPP 值 和 降 
水 .温度 ,太阳 辐射 .实际 蒸 散 量 潜在 蒸 散 量 等 数 
据 , 利 用 主 成 分 分 析 PCA 和 Pearson 相关 系数 法 研究 
了 NPP 变化 的 气候 主 控 因 子 。 

研究 结果 表明 (图 8) ,不 同年 份 不 同 气候 因子 
的 变化 对 研究 区 NPP 贡献 程度 不 一 样 ,但 是 温度 、 
太阳 辐射 .降水 均 是 影响 NPP 变 化 最 相关 的 影响 
子 。 图 8a~ 图 8e 分 别 代表 5 个 不 同时 期 不 同 气 候 因 
子 年 均值 对 NPP 变化 影响 程度 ,1994 一 2018 年 间 温 
度 对 NPP 变化 影响 最 大 ,2000 一 2006 年间 降水 对 
NPP 变化 影响 最 大 ,2012 年 则 是 太阳 辐射 因子 对 
NPP 变化 影响 最 大 。 其 中 ,降水 与 NPP 的 Pearson 相 
关系 数 为 正 值 即 正 相 关 , 温 度 、 太 阳 辐 射 与 NPP 的 
Pearson 相关 系数 为 负 值 即 负 相关 (图 8f)。 由 此 可 
知 , 水 分 是 影响 植物 光合 作用 的 重要 因素 ,降水 量 
的 变化 改变 了 土壤 水 的 含量 ,在 雨水 充足 且 温 度 适 
宜 的 年 份 该 地 区 NPP 显 著 增加 ;由 于 哈萨克 斯 坦 首 


(b) 2018 年 NPP 
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(d) 气候 变化 作用 下 NPP 年 际 变化 
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图 7 气候 变化 作用 下 首都 疾 NPP 变 化 


Fig.7 Changes in NPP in the metropolitan region as a result of climate change 
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注 : 图 f 中 红 框 表示 NPP 和 对 应 
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表示 太阳 辐射 ;temp 表示 平均 温度 。 
图 8 不 同 气 候 因子 主 成 分 分 析 与 Pearson 系数 计算 


Fig. 8 Principal component analysis and determination of the Pearson coefficient for several climatic parameters 


都 圈 位 于 干旱 - 半 干 旱地 区 , 缺 水 情况 下 温度 升 高 
和 大 阳 辐 射 的 增加 导致 水 分 快速 流失 ,从 而 使 植物 
光合 作用 等 生命 活动 产生 了 不 利 影响 。 


4 讨论 


哈萨克 斯 坦 首 都 阿 斯 塔 纳 位 于 干旱 - 半 王 旱 草 
原 区 ,冬季 寒冷 ,夏季 潮 热 ,常年 十 旱 和 大 风 给 当地 
生态 稳定 性 带 来 严重 的 威胁 ,哈萨克 斯 坦 政府 为 改 
善 阿 斯 塔 纳 城市 人 居 环 境 从 1997 年 实施 了 首都 圈 
绿 环 工程 。 为 了 更 好 地 探究 首都 圈 植 被 NPP 的 时 
空 变化 特征 以 及 迁都 对 其 的 影响 ,本 研究 选取 
1994 年 作为 迁都 前 的 对 照 ,每 隔 6 a 提取 一 幅 研 究 
区 土地 利用 覆盖 图 。 研 究 区 范围 的 划 定 依据 为 哈 
萨克斯 坦 提供 的 调查 监测 网 的 。 本 文 在 基于 CASA 
模型 的 基础 上 ,通过 设 定 不 同情 境 , 分 别 计算 研究 
区 内 植被 NPP ,探究 在 人 类 活动 和 气象 条 件 变 化 下 
对 研究 区 内 植被 NPP 的 影响 ,最 后 通过 Pearson 相关 
系数 计算 和 PCA 分 析 确 定 研 究 区 内 植被 NPP 变 化 


的 驱动 因素 。 

人 类 活动 是 中 亚 地 区 植被 初级 生产 力 动态 变 
化 的 重要 驱动 力 ™ ,研究 发 现 从 1994 一 2018 年 哈 萨 
克 斯 坦 首 都 圈 绿 环 工程 的 实施 使 大 面积 的 草原 和 
其 他 植被 转换 为 森林 ,植被 NPP 整体 呈现 上 升 的 趋 
势 ( 图 3)。1994 一 2000 年 间 哈 萨克斯 坦 将 首都 迁 至 
于 此 , 短 时 间 内 人 口 的 增长 建筑 面积 扩 增 、 用 水 消 
耗 过 度 导 致 地 表 水 分 减少 ,水 分 减少 降低 了 草原 植 
物 的 光合 作用 效率 ,抑制 植物 活动 和 有 机 质 生产 ， 
最 终 降 低 了 植被 生产 力 ”。 随 着 “ 绿 环 工程 ”的 实 
施 和 后 续 地 表 水 逐步 恢复 ,森林 覆盖 面积 不 断 扩 
大 ,研究 区 内 NPP 也 在 不 断 的 增加 。 男 外 ,在 2000 
年 以 后 ,哈萨克 斯 坦 经 过 相应 的 政策 改革 ,开始 恢 
复 农 牧 业 ,使 得 该 地 区 生产 规范 化 2 , Han SE 
究 表明 干旱 地 区 合理 的 放牧 能 有 效 防止 植被 退化 ， 
有 利于 植被 生态 系统 的 可 持续 发 展 。 

半 干 旱 - 干 旱 生 态 系 统 中 的 植被 对 气候 变化 极 
其 敏感 ,降水 是 影响 地 区 植被 NPP 的 主要 因素 。 黄 


202301.00134v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1926 干 F Wt F 39 卷 
钰 “发 现 ,中国 陆地 植被 NPP 变 化 主要 受降 水 量 的 。 更 强 。 


主导 作用 ; 韩 其 飞 等 呈 的 研究 也 表明 ,干旱 地 区 降 
水 量 是 影响 NPP 变 化 的 主要 因子 ;而 植被 NPP 变 化 
是 基于 多 种 因素 共同 影响 的 结果 ”*””W。 本 文 为 量化 
人 类 活动 和 气候 变化 对 地 区 植被 NPP 产 生 的 影响 ， 
模拟 了 气候 条 件 保 持 在 1994 年 水 平 不 变 的 情境 下 ， 
计算 了 不 同年 份 植被 NPP( 图 5)。 结 合 图 3d 可 以 看 
出 ,在 人 类 活动 的 影响 下 研究 区 内 植被 NPP 增 长 速 
率 最 快 。 其 次 模拟 了 土地 覆盖 和 NDVI 值 保持 为 
1994 年 的 水 平 不 变 ( 图 7), 可 以 看 出 气候 变化 对 人 研 
究 区 内 植被 NPP 影 响 波动 较 大 ,不 同年 份 之 间 NPP 
差 值 差异 明显 。 研 究 表 明 : 降 水、 太阳 辐射 和 温度 
是 哈萨克 斯 坦 首都 圈 植 被 NPP 变化 的 主要 驱动 因 
R ,水 分 是 影响 植物 光合 作用 的 重要 因素 。 降 水 量 
的 变化 改变 了 土壤 水 的 含量 ,因此 在 雨水 充足 且 温 
度 适 宜 的 年 份 NPP 显 著 增 加 ,水 热 良 好 搭配 在 NPP 
形成 中 起 到 重要 的 作用 YW。 但 是 由 于 哈萨克 斯 坦 
首都 圈 处 于 干旱 -半生 旱地 区 ,在 缺 水 的 情况 下 , 温 
度 升 高 和 太阳 辐射 增加 会 导致 水 分 快速 流失 及 植 
被 气孔 关闭 ,同时 伴随 着 潜在 燕 散 和 水 分 胁迫 增 
强 , 进 而 导致 草地 生产 力 降 低 ,从 而 对 植物 光合 作 
用 等 生命 活动 产生 不 利 影响 *。 


5 结论 


本 文通 过 CASA 模 型 模拟 了 3 种 不 同情 境 下 哈 
萨克斯 坦 首都 轿 植 被 净 生 产 力 的 时 空 变化 ,定量 分 
析 了 人 类 活动 和 和 气候 变化 对 该 地 区 NPP 的 影响 , 探 
讨 了 不 同 气候 因子 对 NPP 的 驱动 力 大 小 , 主要 结论 
如 下 : 

(1) 在 1994 一 2018 年 间 ,哈萨克 斯 坦 首都 圈 植 
被 NPP 呈现 波动 上 升 趋势 ,NPP 多 年 平均 值 为 
226.21 g Cm >”a'。 首 都 轿 5 种 主要 十 地 覆盖 类 型 
NPP 均 值 从 高 到 低 分 别 为 : 防护 林 > 其 他 值 被 > 裸 
地 > 建筑 用 地 > 水 体 。 

(2) 人 类 活动 对 首都 圈 NPP 影 响 有 显著 增益 效 
果 (0.38 Tg C- a`, P<0.01) ,气候 变化 对 首都 圈 NPP 整 
体 上 呈 负 增益 效果 (-0.07 Tg Ca”, P=0.34), AŽ 
活动 (土地 覆盖 变化 ) 和 气候 变化 共同 作用 下 PB 
克 斯 坦 首 都 圈 NPP 总 体 上 呈现 显著 正 增 益 效 果 
(0.27 Tg Ca ,P<0.1)。 相 较 于 气候 变化 ,人 类 活动 
对 植被 净 初 级 生产 力 影 响 更 加 明显 ,相关 性 也 


(3) 气候 因子 中 降水 温度 ,太阳 辐射 和 实际 蒸 
散 量 是 影响 研究 区 NPP 变 化 的 直接 驱动 因素 ,其 中 
降水 对 研究 区 NPP 有 显著 的 正 效应 ,温度 ,太阳 辆 
射 和 实际 蒸 散 量 对 研究 区 NPP 有 显著 负 效 应 。 
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Characteristics and drivers of the spatial-temporal change of net primary 
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Abstract: Clarifying the impacts of climate change and human activity on plant’ s net primary productivity 
(NPP) is crucial for regional ecosystem transformation and sustainable development. Using a simulated scenario 
experimental design based on the CASA model, five periods of Landsat remote sensing images, and 
meteorological data spanning the years 1994 to 2018, this paper determines the effects of two factors, human 
activities and climate change, on plant’ s NPP in the capital area. The results indicate that: (1) During 1994-2018, 
the multi-year average of NPP in the studied area was 226.21 g C-m”°-a`', with a fluctuating rising trend; (2) The 
execution of the Green Ring Project led to a large positive gain in NPP owing to an increased plantation area 
(0.38 Tg C:a ', P < 0.01).In contrast, the impact of climate change on NPP is more variable, with an overall loss 
effect (-0.07 Tg C-a', P = 0.34). Under the combined influence of human activity (land use change) and climate 
change, the NPP of Kazakhstan’ s capital area exhibited a considerable positive gain effect (0.27 Tg Ca ', P < 
0.1); climate change has a lesser impact on plants’ net primary productivity than human activities. (3) 
Temperature, solar radiation, and precipitation are the most influential climatic elements on NPP. During 1994- 
2000 and 2006-2012, increased temperature and reduced precipitation caused NPP loss, which dropped from 
218.50 g Cem” -+a ', 201.19 g C-m”-a' to 189.00 g C-m a’', 188.48 g C-m a; With the improvement of 
precipitation circumstances over 2000-2006 and 2012-2018, the mean value of NPP in this area increased 
significantly, reaching 201.19 g C-m~*+a ', 207.73 g C:m av correspondingly. These research findings assist in 
elucidating the processes of climate change and human activities on NPP and may also serve as a guide to 
enhance the ecological quality of the desert-steppe area, alleviate the global warming issues, and serve the carbon 
neutrality goal of Kazakhstan by 2060. 

Keywords: CASA model; NPP; Kazakhstan capital circle; PCA; simulation experiment design 


